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Axial  strain  localization  in  a  mesoscopic  chain  of  five  bistable 
springs.  The  spring  collapse  mimics  the  local  kinking  of  com¬ 
pressed  carbon  nanotubes. 
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Predicted  stress-strain  response  (solid  line)  at  the  macroscopic 
scale,  reproducing  the  experimental  behavior  of  a  real  CNT 
foam  (dashed  line). 
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•  Axial  strain  localization  in  microscopic  bistable  spring  chains  mimics  kink¬ 
ing  of  compressed  carbon  nanotube  arrays.  •  Infinitesimal  viscous  events  at 
the  microscale  induce  time-independent  hysteresis  at  the  mesoscale.  •  Multi- 
scale  mechanical  modeling  of  CNT  foams  is  obtained  through  an  information- 
passing  approach.  •  Available  experimental  results  on  compressed  CNT 
foams  are  reproduced  with  excellent  agreement. 
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scale  model  formed  by  a  chain  of  nonuniform  mesoscopic  springs  is  employed 
to  describe  the  entire  CNT  foam.  A  numerical  approximation  illustrates  the 
main  features  of  the  proposed  multiscale  approach.  Available  experimental 
results  on  the  compressive  response  of  CNT  foams  are  fitted  with  excellent 
agreement. 
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1.  Introduction 


Since  their  discovery  ( Radushkevich  and  Lukyanovich 

BUM 


1952  Oberlin  et 


Iijima  et  ah,  1995),  carbon  nanotubes  (CNTs)  have  been  widely 


studied  for  the  understanding  of  different  aspects  of  their  chemical,  electri¬ 
cal  and  mechanical  responses.  Because  of  their  unique  properties  and  their 
multiscale  nature,  forests  of  vertically  aligned  carpets  of  CNTs  have  been 


proposed  for  a  variety  of  applications  (Gjerde  et  al.  1991;  Veedu  et  al. 


Daraio  et  al. 

2004b; 

Maj umder  et  al. 

2005 

Maheshwari  and  Saraf, 

2008) 

The  mechanical  response  of  individual  nanotubes  under  axial  and  rac 


2006 


ial  de¬ 


formation,  and  their  bending/buckling  modes,  have  been  studied  extensively 


using  experimental,  theoretical  and  molecular-dynamics  analysis  (Ii 

jima  et 

al.,  1995|;  Yakobson  et  al. 

1996 

Falvo  et  ah,  1997|  Belytschko  et  al 

,  2002; 

Arroyo  and  Belytschko 

2003)  I 

-’antano  et  al. 

2004 

Cao  and  Chen 

2006). 

The  elastic  modulus  E  of  individual  carbon  nanotubes  has  been  reported 
to  be  very  high:  ~  1  TPa  (Pantano  et  ah,  2004).  However  in  experiments 


this  value  can  vary  widely,  depending  on  the  number  of  defects,  the  CNT 
microstructure  and  the  synthesis  method  followed. 

The  study  of  the  mechanical  properties  of  CNTs  was  later  extended  to 


bundles  of  nanotubes  under  pressure  (C 

Uiesnokov  et  ah, 

1999 

Peters  et  al. 

2000  Chan  et  al.  2003 

Liu  et  al. ,  2005 

Qi  et  al. 

2003) 

and  to  CNT  forests 

under  nanoindentation 

Mesarovic  et  al. 

2007 

.  The  study  of  thin  structural 

foams  (Gibson  and  Ashby, 

1998)  for  cushioning  ( 

Zhang  et  al. 

2009),  energy 

dissipation  (Teo  et  ah, 

200 

7)  and  protection 

Liu  et  al. 

2008 

)  has  recently 

received  increasing  attention  for  several  practical  applications,  including  mit¬ 
igation  of  explosive  loading  (Nesterenko,  2001).  Nanotube  based  films  have 


been  reported  as  an  excellent  alternative  to  regular  foams,  exhibiting  a  super- 


compressible  foam-like  behavior  under  compressive  cycling  loads  (Suhr  et  al. 


2007 

Teo  et  al. 

2007 

Tao  et  ah, 

2008 

Deck  et  ah, 

2007 

Cao  et  al. 

2005 

2 


Investigations  on  the  dynamic  response  of  foam-like  forests  of  CNTs  under 


dynamic  ball  impacts  have  also  been  performed  ( 

Daraio  et  al. 

2004b 

Daraio 

et  al. 

2004b; 

Daraio  et  al. 

2006 

Misra  et  al. , 

2009 

) .  Results  show  a  strongly 

nonlinear  response  that  appears  to  be  very  suitable  for  energy-absorbing  lay¬ 
ered  materials  in  noise  and  shock  wave  mitigation  and  as  nonlinear  springs 
for  assembling  nonlinear  acoustic  crystals.  The  CNT  forests’  response  was 
also  found  to  be  strongly  dependent  on  the  forests’  microstructure  (height, 
density,  alignment,  etc.)  and  growth  method.  In  certain  cases  the  possi¬ 
ble  presence  of  plastic  deformation  and  fracturing  of  the  tubes  was  reported 


(Daraio  et  al.  2004b). 


Mechanical  models  consisting  of  chains  of  identical  bi-stable  springs  have 
been  extensively  studied  by  several  authors,  since  the  can  describe  a  se¬ 
ries  of  relevant  nonlinear  material  behaviors  (e.g.,  phase  transformations, 
reversible  pseudo-plasticity,  hysteresis,  fracture),  through  the  interplay  be¬ 
tween  macroscopic  and  microscopic  length  scales  (refer  e.g.  to  Ortiz,  [1999 


Puglisi  and  Truskinovsky  2000  2002  [2005).  It  is  well  known  that  such  sys¬ 


tems  exhibit  a  ‘bumpy’  multi-well  energy  landscape,  allowing  for  multiple 


metastable  equilibria  (cf.,  e.g., 

Blesgen, 

2003 

Braides  and  Cicalese 

2007[ 

Charlotte  and  Truskinovsky,  2002,  2008 

Ortizi  19991  Puglisi  and  Truski- 

novsky 

2000, 

2002).  In  particular, 

Pampolini  and  Del  Piero 

(2008 

)  have 

recently  found  that  they  well  describe  the  hysteretic  response  of  open-cell 
polyurethane  foams  under  confined  compression  tests. 

In  this  article,  we  present  a  phenomenological  model  of  the  mechani¬ 
cal  response  of  carbon  nanotube  foams  under  compressive  loading,  which 
is  inspired  by  some  distinctive  features  of  the  micromechanical  response  re¬ 


lumped  masses  connected  by  nonlinear  springs.  It  captures  the  ‘three-phase’ 
compressive  deformation  response  of  CNT  forests  shown  by  a  number  of  ex¬ 
perimental  studies.  The  compressive  deformation  response  is  characterized 
by  an  initial  elastic  deformation,  a  non-homogeneous  buckling  (or  plateau) 
regime,  often  featuring  a  sawtooth-like  profile,  and  a  densiheation  phase. 


This  three-phase  response,  which  is  common  for  cellular  materials  (Gibson 


and  Ashby,  1998),  is  usually  accompanied  by  marked  hysteresis  and  strain 
localization  in  CNT  structures.  We  show  in  Section  |3]  that  a  series  of  bistable 
elastic  springs  (Fig.  [l|  described  by  the  potential  in  Eq.  ([TJ)  exhibits  a  simi¬ 
lar  stress-strain  response  (Fig.  [6]),  and  through-the-thickness  localization  of 
the  axial  deformation  (Fig.  [7]).  The  latter  effect  mimics  the  snap-buckling 


3 


events  observed  through  Scanning  Electron  Microscope  (SEM)  in  real  CNT 
arrays  (cf.  Fig.  2  of  Cao  et  al.  2005;  Fig.  3  of  Hutchens  et  al.|  2010).  Such 


a  model  therefore  appears  to  be  effective  in  describing  the  microstructure 
rearrangements  taking  place  in  compressed  CNT  foams,  within  a  simple  ID 
framework. 

We  focus  on  the  time-independent  component  of  the  hysteresis  associated 
with  the  compressive  loading/unloading  of  CNT  foams,  which  is  essentially 
due  to  effects  such  as  friction,  entanglement  and  electrostatic  interaction  be¬ 
tween  individual  and  bundles  of  carbon  nanotubes  (Suhr  et  ah,  2007; 


et  al. 

2007 

al., 

20( 

15). 

Teo 


Tao  et  al.  2008  |Deck  et  al.  2007  Misra  et  al.  2009|  Cao  et 


al.[  2005).  The  macroscopic  hysteresis  is  described  as  a  rate-independent 
phenomenon  induced  by  the  succession  of  infinitesimal  viscous  events  at  the 
microscopic  scale  (Puglisi  and  Truskinovsky ,  2005).  As  a  result,  the  me¬ 


chanical  model  presented  in  this  work  does  not  account  for  viscosity  or  other 
rate-dependent  effects  at  the  macroscopic  scale,  which  we  address  to  future 
work. 

Most  of  the  available  studies  on  bistable  spring  chains  consider  only  two 
spatial-temporal  micro-macro  scales.  The  formulation  adopted  in  the  present 
study  instead  introduces  three  different  time-space  scales:  a  ‘microscopic’ 
scale  (of  the  order  of  nanometers),  which  is  associated  with  the  individual 
bi-stable  springs  and  an  infinitesimal  portion  of  the  total  foam  thickness  Ltot ; 
a  ‘mesoscopic  scale’  (of  the  order  of  micrometers),  corresponding  to  the  limit 
of  an  infinite  series  of  microscopic  springs  and  representing  a  finite  portion 
of  Ltot ;  and  a  ‘macroscopic’  scale  (of  the  order  of  millimeters),  describing 
the  entire  foam.  We  start  with  the  derivation  of  a  discrete  ID  model  at 
the  microscopic  level,  showing  N  +  1  particles  having  nearest  neighbors  con¬ 
nected  by  N  uniform  nonlinear  springs.  The  discrete  spring  potentials  are 
chosen  to  allow  the  study  of  softening  and  strain  localization  of  the  tubes. 
This  is  discussed  in  Section  [2}  Here,  we  build  on  the  concept  of  bi-stable 


springs  discussed  by  Puglisi  and  Truskinovsky  in  ( 

Puglisi  and  Truskinovsky| 

2000  Puglisi  and  Truskinovsky,  2002 

Puglisi  and  Truskinovsky 

2005  ).  In 

Section  [3]  we  determine  analytically  the  continuum  limit  N  — »  oo  of  the  mi¬ 
croscopic  spring  chain,  by  particularizing  to  the  present  case  the  analysis 
presented  in  Puglisi  and  Truskinovsky  (2005).  Such  a  limit  defines  a  meso¬ 
scopic  dissipative  element.  In  Section  [4]  we  then  formulate  an  upper  scale 
model  of  the  entire  CNT  foam  through  an  information-passing  approach, 
by  superimposing  a  finite  number  of  mesoscopic  dissipative  springs  with  dif¬ 
ferent  mechanical  properties.  We  account  in  this  way  for  inhomogeneities 


4 


induced  by  the  CNT  growth  process,  which  several  authors  think  be  a  lead¬ 


ing  cause  of  the  discrete  folding/buckling  events  described  above  (Cao  et  al. 


2005  Hutchens  et  al.  2010).  The  introduction  of  non-uniform  mechanical 


properties  allows  us  to  enrich  the  formulation  given  in  Puglisi  and  Truski- 


novsky 

(2000 

), 

Puglisi  and  Truskinovsky 

(2002) 

Puglisi  and  Truskinovsky 

(2005), 

Pampolini  and  Del  Piero 

, (2008 

)  for  bistable  spring  chains,  modeling 

macroscopic  hardening,  instead  of  a  perfectly  ‘plastic’  response.  Hardening- 
type  post-buckling  regimes  are  commonly  observed  in  compression  tests  on 
CNT  foams  (Cao  et  al.  2005  Misra  et  al.  2009).  We  present  in  Section  [5] 


a  numerical  micro-meso  convergence  study,  and  the  fitting  of  mesoscopic 
model  parameters  to  different  available  experimental  data  on  compressed 
CNT  foams.  We  demonstrate  that  the  proposed  model  and  its  numerical  im¬ 
plementation  are  capable  of  reproducing  the  physical  behavior  of  real  CNT 
foams  with  excellent  agreement.  We  end  the  article  with  a  critical  evaluation 
and  discussion  of  the  results  and  an  outlook. 


2.  The  mechanical  model  at  the  microscopic  scale 

We  model  an  infinitesimal  portion  of  a  CNT  array  as  a  collection  of  N  + 1 
lumped  masses  m°, . . .  ,mN  piled  up  one  over  the  other  (with  N  >  2).  In 
this  configuration,  m°  is  clamped  at  the  bottom  of  the  pile  at  the  position 
£°  =  0,  whereas  mN  is  on  top  at  position  xN  =  L  >  0,  referring  to  the 
unstressed  reference  configuration  (Fig.  [lj^).  The  nearest  neighboring  mass 
points  are  connected  by  N  nonlinear  microscopic  springs.  We  assume  that 
the  reference  configuration  displays  equal  distances  Hn  L/N  between  the 
masses,  and  denote  the  axial  displacement  of  the  mass  ml  (positive  upward) 
by  ulN  (with  u°N  =  0,  see  Fig.  [jj).  We  set  u ^  :=  {u°N, 

For  the  mechanical  energy  V1  of  the  microscopic  spring  placed  between 
nearest  neighbors  m!_1  and  m!,  we  assume  the  three-branch  expression  de¬ 
fined  by 


4(4) 


'  4(4)  =  -k ;(,[4  +ln(l-e<)],  4  <  4, 

<  4(4)  =  ci  +  44  + -44  4  <  4  <  4> 

k  4(4)  =  c2-4[4-4+in(i-(4-4))]>  4  <  4, 


(i) 
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Figure  1:  (Color  online)  (a)  Schematic  diagram  of  a  vertically  aligned  CNTs  foam,  uni¬ 
formly  loaded  in  compression,  (b)  SEM  of  the  as  grown  carbon  nanotube  film  showing  the 


alignment  and  the  microstructural  layering  due  to  the  growth  process  ( Deck  and  Vecchio 


2005 ) .  (c)  Modeling  of  a  portion  of  a  CNT  foam  as  a  collection  of  microscopic  mass-spring 


elements. 


with  £l  =  £l  —  e q,  where 


£l  =  £\uN ) 


U 


.i—  1 
N 


U 


N 


h 


N 


(2) 


is  the  strain  measure  associated  with  such  a  spring  (positive  in  compression). 
In  ([!]),  the  quantity  £l0  >  0  determines  the  value  of  £l  corresponding  to 
the  first  minimum  of  V1  (‘equilibrium’  or  ‘initial’  strain);  kl0  >  0,  klb  <  0, 
klc  >  0,  £la  >  0  and  £lc  >  £la  are  constitutive  parameters  (five  independent 
parameters);  the  constants  Ci  <  0  and  c2  >  0  are  such  that  V*(£la)  =  (£a)> 

Vbz(il)  =  Vc(il);  and  it  results 
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a: 


el  =  el 


ki  +  <7®  ’ 


(3) 


with 


4 


4  (4  + 


ki  + 


crl 


+ 


a:  —  a: 


)  {ki  +  ejkj  +  ejaj) 


(K  +  oj)  (ki  +  ai) 


(4) 


<4  =  fcn 


1  —  £l 


—  aa  +  H\£c~  £a)- 


(5) 


The  relationship 


ut 

i_ei  > 


£*  <  el. 


<y{(4)  =  v1'  =  {  4 +  4(4 -4),  4<4<e 


o 


kU^-ej) 
k  !-(£*-£*)’ 


4  <4 


(6) 


defines  the  stress  a1  acting  in  the  spring  connecting  ml  1  with  rri 
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We  deduce  from  ([I])  that  V1  owes  the  two- well  prohle  shown  in  Fig.  [2j  Ac¬ 
cordingly,  the  stress-stain  relationship  (J6|  is  described  by  the  non-monotone 
prohle  depicted  in  Fig.  [3j  It  is  worth  noting  that  the  mechanical  response 
of  the  generic  microscopic  spring  encompasses  two  stable  phases,  for  £*  <  e\ 
(phase  a)  and  e1  >  (phase  c),  respectively;  and  an  intermediate  unstable 
phase  b  for  ela  >  (spinodal  regime).  For  future  use,  we  set 


An*  :=  d  -  ai 


e\ 


£  <J 

n'-'  n 


ai  +  ei  Act* 


(7) 


and  let  alM  denote  the  slope  of  the  linear  branch  of  the  convex  hull  of  V1 
(Maxwell  stress,  see  Fig.  [3]).  We  also  convene  to  denote  the  value  of  el  where 
a1  =  alM  in  phase  a  by  el  ,  and  the  value  of  e1  such  that  a1  =  alM  in  phase  c 
by  £*m.  With  reference  to  the  generic  spring,  Equation  (J6|i  highlights  that 
klQ  represents  the  elastic  stiffness  at  zero  stress  in  phase  a  (s*  =  0).  On  the 
other  hand,  Equation  points  that  klb  represents  the  (constant)  stiffness 
in  the  spinodal  phase  b ,  while  Equation  ([6J3  reveals  that  klc  represents  the 
stiffness  at  zero  stress  in  phase  c  (e*  =  £*).  We  name  ‘symmetric’  the  case 
with  klc  =  kl0  and  ‘asymmetric’  the  one  with  fc*  7^  kl0.  The  meaning  of  the 
other  quantities  appearing  in  Equations  (Jl])  -  ([7])  is  illustrated  in  Figs.  [2]  and 

IS 


3.  Dynamic  relaxation  and  hysteresis  at  the  mesoscale 

There  are  two  fundamental  characteristics  reported  experimentally  in  the 
bulk  compressive  response  of  dense  vertically  aligned  CNTs  forests:  (i)  a 
good  recovery  of  deformation  even  at  large  compressive  strains  after  a  suf¬ 


ficiently  large  recovery  time  (Cao  et  al.  2005  Misra  et  al.  2009);  and 
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a  strong  hysteresis,  observed  in  particular  in  as-grown  foams,  attributed  to 
several  potential  effects,  including  friction,  entanglement  and  electrostatic 
interaction  between  individual  and  bundles  of  carbon  nanotubes  (Suhr  et  ah, 


2007 

Teo  et  al. 

2007 

Tao  et  al. 

2008 

Deck  et  al.  2007;  Misra  et  al. 

,  2009 

Cao  et  al. 

2005 

).  In  the  present  Section  we  will  analyze  such  a  dissipative 

behavior,  by  studying  a  dynamic  switching  process  at  the  microscopic  scale 
between  the  phases  (a)  and  (c)  described  in  Figs.  [51  and  |3j  This  is  in  line 


with  the  ideas  in  1 

3uglisi  and  Truskinovsky 

(5 

>002). 

Puglisi  and  Truskinovsky 

(2005 

) .  Following 

Puglisi  and  Truskinovsky 

(200£ 

i),  we  name  a  response  of 

the  material  plastic ,  if  the  strain  el  of  a  single  spring  exceeds  el0  +  ela.  For 


V1 


Figure  2:  (Color  online)  Mechanical  energy  V1  of  the  generic  microscopic  (bi-stable)  spring. 


a  chain  of  N  springs,  this  can  be  characterized  by  the  occurrence  of  loading 
and  unloading  stress  plateaux.  However,  we  notice  that  our  analysis  excludes 
accumulation  of  permanent  deformation,  therefore  the  end  point  of  one  hys¬ 
teresis  cycle  coincides  with  the  start  point  of  the  next  cycle.  In  this  sense,  the 


present  hysteresis  model  is  time-independent.  We  will  show  in  Section  5.2 


that  it  is  nevertheless  capable  of  capturing  a  previously  accumulated  perma¬ 
nent  deformation  ( mechanical  preconditioning ),  through  suitable  definition 
of  the  initial  strains  elQ. 

Within  the  current  Section,  we  rescale  for  simplicity  L  to  unity,  and 
assume  that  V1  is  independent  of  the  spatial  position.  Accordingly,  we  drop 
the  superscript  i  in  front  of  the  spring  properties.  This  because  we  refer 
the  following  analysis  to  a  finite  portion  of  the  CNT  foam,  regarding  such 
a  mesoscopic  element  as  the  limit  for  N  — »  oo  of  a  series  of  N  identical 
microscopic  springs.  Furthermore,  we  restrict  our  attention  to  the  case  with 
kc  =  ko  (’’symmetric”  model),  and  Eq  =  0.  For  the  present  analysis  we 
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Figure  3:  (Color  online)  Stress  a 1  versus  strain  el  relationship  in  the  generic  microscopic 
spring. 


require  a  certain  smallness  condition  on  ea  and  ec  relating  to  strong  pinning 


that  disappears  in  the  limit  N  — »  oo,  see  Puglisi  and  Truskinovsky  (2002) 


We  define  the  mechanical  energy  of  the  foam  as 


En{un )  -  —  ^\d(e*(MAr))  (8) 

i= 1 

with  the  effective  potential 

{-/c0[£  +  ln(l  -  £)]  if  £  <  ea, 

ci  +  aae  +  -  Ea)2  if  £a<e<£C) 

c2  -  k0[e  -  e*  +  ln(l  -  c  +  e*)]  if  ec  <  e. 

Let  a  be  the  given  total  stress,  coinciding  at  equilibrium  with  the  stress  in 
each  individual  spring  (a  =  a1  =  ...  =  crN).  The  total  average  strain  is 
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simply 


1  N 

e{uN)  :=  —  ^V(V). 

i=l 


We  model  plasticity  by  the  gradient  flow  equations  (Puglisi  and  Truski- 
vi*(uN)  = 


novsky  2005) 


,  i 


(e^ujv), . .  .,eN(uN)) 


dei 


(9) 


with  the  total  energy 


1  N 

i= 1 

The  evolution  equation  (|9])  lets  el  evolve  towards  local  minimizers  of  $at. 
We  are  interested  in  the  limit  v  —*  0  which  amounts  to  infinitely  fast  evolution 
such  that  s(u]sr)  attains  a  local  minimizer  of  Q]y.  First  we  construct  the 
equilibrium  points.  Inside  the  i-th  spring  element,  the  strain  must  satisfy 
the  condition  V\el)  =  cr.  For  given  total  stress  cr,  there  are  at  most  the 
three  local  minimizers 

^  ^  _  cr  —  <r(l  +  e*)  +  k0e*  _  „ 

~r  |  ;  T  &ai  T  £*• 

k0  +  cr  kb  ko  +  cr 


Let  p,  q,  1  —  p  —  q  denote  the  phase  fractions  of  the  minimizers  a,  b,  and  c, 
which  corresponds  to  having  Np,  Nq,  IV (1  —  p  —  q)  springs  in  phase  a,  b,  and 
c,  respectively. 

As  £  t— >  V (e)  is  concave  in  Regime  b,  if  the  elongation  of  a  spring  in  the 
local  minimum  £b  is  altered  by  an  arbitrarily  small  perturbation,  it  will  move 
(according  to  the  sign  of  the  perturbation)  to  either  £a  or  ec.  In  consequence, 
any  system  of  N  springs  with  q  ^  0  is  unstable  and  we  may  in  the  following 
restrict  to  the  case  q  —  0. 

From  £  =  p£a  +  (  1  —p)£c  we  compute  the  equilibrium  stress-strain  relation 
to  be 


a(£) 


fc0(g  -  ep) 

1  -  (e  -  £P) 


(10) 


with 


£p  :=  (1  -p)g* 


that  can  in  a  natural  way  be  identified  with  the  plastic  strain.  From  (10)  we 
see  that  cr  only  depends  on  the  elastic  strain  gei  :=  £ 


-p- 
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For  the  energy  of  the  equilibrium  configuration  with  Np  springs  in  phase  a 
we  find 


Ep(e)  =  —h 


a(e) 


+  In 


+  (1  -p)c2 


. k0  +  cr(£ )  Vfc0  +  er(e)  /  . 

-k0  [(e  -  ep)  +  ln(l  -  (e  -  ep))]  +  (1  -  p)c2. 


Note  that  Ep  has  a  finite  number  of  local  minimizers  (depending  on  the 
remaining  parameter  p).  As  explained  in  Puglisi  and  Truskinovsky  (2005), 
the  switching  takes  place  between  branches  that  differ  in  exactly  one  element 
in  the  phase  state  and  the  succession  of  N  such  steps  describes  the  transition 
from  one  homogeneous  state  to  the  next.  Each  of  these  steps  can  be  thought 
of  as  the  combination  of  an  elastic  part  and  a  plastic  part. 

The  stress-strain  curve  of  the  foam  follows  a  sawtooth  pattern  as  illus¬ 
trated  in  Fig.  [4]  (left).  We  denote  by  A 1  the  end  point  of  the  i-th  branch. 
The  A1  are  the  final  states  of  the  elastic  steps  ( Bl  — >  Al+l)  where  the  system 
remains  on  the  same  metastable  branch  as  long  as  possible.  The  plastic  steps 
(. A 1  — >  Bl)  are  characterized  by  that  the  total  strain  is  fixed  and  the  system 
switches  between  metastable  branches  that  are  neighbors  ([Np]  =  1,  and  [Q\ 
generically  denotes  the  jump  of  a  quantity  Q).  These  considerations  lead  to 
the  representation 


Puglisi  and  Truskinovsky  (2005 


Ai 


fcpgq  \ 

1-J’ 


Bi 


f  k0(£a-  £*/N)\ 


1  <  i  <  N 


in  the  (e:,  cr)-diagram,  where 


i  —  1 

£Ai  -=  £a  H 


(11) 


So  we  can  compute  that  a  plastic  step  is  characterized  by 


[e]  =  0,  [Np]  =  1,  [or] 
whereas  an  elastic  step  fulfills 

[e]  =  %,  m  =  0,  M 


k{)£ * 

N(l-ea)(l  +  e*/N-eay 


k0£* 

N(1  —  ea)(l  +  £*/N  —  £a) 
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a 


Figure  4:  Overall  stress-stain  pattern  for  a  chain  with  TV  =  10  microscopic  springs  (left), 
and  limit  response  for  N  — >  oo  (right). 


The  evolution  equation  (J9])  lets  £l  evolve  towards  local  minimizers  of 
Now  we  want  to  look  at  the  energetics  of  the  plastic  and  the  elastic  regime. 
For  an  elastic  step  we  have  the  energy  difference 

AEn  =  fc0[ln(l  -  £a)  -  ln(l  -  £a  T  £*/iV)]  +  k°£*  +  °2 , 


In  the  same  spirit,  we  calculate  that  for  a  system  with  N  >  1  springs,  the 
plastic  dissipation  is 


Dm  kn 


ln(l  -  £a  +  e*/N)  -  ln(l  -  £a) 


c2  +  k(j£^ 

N 


,  1  £* 
k°^N 


c2  + 
N 


for  a  £n  G  (1  -  £a,  1  -  £a  +  £*/N).  (12) 


Clearly,  £N  —> ►  1  —  £a  for  — >■  oo. 

In  one  hysteresis  cycle,  there  are  N  loading  steps  and  N  steps  when  the 
system  is  unloaded,  so  we  have  totally  2 N  steps  that  dissipate  energy.  The 
total  dissipated  energy  D  in  a  cycle  becomes  in  the  limit  N  — >  oo 


D 


lim  2 NDn  = 

TV— >oo 


‘2h()£^£a 

1  &CL 


-  2  c2 


(13) 
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We  put  the  modulus  here  to  ensure  that  the  dissipation  is  positive. 

The  limit  stress-stain  pattern  for  N  — >  oo  is  shown  in  Fig  |4j  It  corre¬ 
sponds  to  a  ’’perfectly  plastic”  behavior  with  stress  plateaux  at  a  =  oa  = 
k0ea/(  1  —  ea )  (loading  plateau)  and  a  =  ac  =  aa  +  Act  (unloading  plateau). 
We  emphasize  again  that  this  ansatz  only  works  for  rate-independent  plas¬ 
ticity  where  the  energy  only  depends  on  start  point  and  end  point  of  the 
evolution,  but  not  on  the  evolution  path  itself.  The  limit  dissipation  (13) 
equals  the  area  enclosed  by  the  limit  stress-strain  response.  As  already  ob¬ 
served,  we  confer  the  behavior  shown  in  Fig.  [4]  (right)  to  a  mesoscopic  spring 
element,  which  represents  a  finite  portion  of  the  CNT  foam  thickness. 


4.  Multiscale  numerical  modeling 

We  formulate  in  the  present  Section  a  multiscale  numerical  model  of  a 
nonlinear  mass-spring  chain,  where  each  spring  represents  either  a  micro¬ 
scopic  bistable  element  (cf.  Section  [2]),  or  a  mesoscopic  dissipative  element 
of  the  kind  introduced  in  the  previous  Section.  We  introduce  two  different 
time  scales:  an  external  (slow)  time  r  G  [to,  t\\  ruling  an  evolution  law  of  the 
applied  boundary  conditions,  and  an  internal  (fast)  time  t  G  [ioAi]  governing 
the  dynamic  relaxation  of  the  system  for  fixed  r.  Depending  on  the  adopted 
model  for  the  individual  springs,  we  may  have  the  t  corresponds  to  the  mi¬ 
croscopic  timescale  (i.e.  to  the  time  ruling  the  microscopic  behavior)  and  r 
to  the  mesoscopic  time  (micro- meso  transition),  or,  alternatively,  that  t  cor¬ 
responds  to  the  mesocopic  time  and  r  to  the  macroscopic  time  (meso-macro 
transition). 

Let  us  denote  a  prescribed  displacement  time-history  of  the  topmost  mass 
mN  by  S(t).  We  introduce  a  discretization  {ti,  . . . ,  tm}  of  the  loading  interval 
[to,ti]  and  compute  the  system  response  for  fixed  r  =  77  and  5  =  5(rk), 
through  integration  with  respect  to  t  of  the  evolution  equations 


—  4+1 


m‘  uN  +  7'  uN  —  a 


-  <7 


i  =  1, 


,N, 


(14) 


which  generalize  the  gradient  flow  equations  (J9])  (Puglisi  and  Truskinovskyj 
2005).  In  (14),  ulN  =  ulN(t )  denote  transient  displacement  histories  of  the 
masses  m°, . . . ,  mN  at  fixed  r;  a1  indicates  the  current  stress  in  the  ith  spring 
(it  is  understood  that  it  results  a1  =  0  for  i  >  A);  and  71,. . .  ,7^  denote 
damping  coefficients. 

Since  we  are  only  interested  in  the  final  equilibrium  configuration  of  the 
transient  internal  motion,  we  ‘overdamp’  such  a  motion  down,  by  introducing 
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fictitious  masses  and  damping  factors  in  (14).  In  detail,  we  set  the  integration 


time  step  At  to  unity,  and  introduce  fictitious  masses  rnl  =  akl  At2,  with 
ki  =  +  Vl+1)"  and  a  >  100.  This  ensures  At  <  0.1  \Jmi /kl  (i  = 

1, ...,  N)  (Fraternali  et  ah,  2009).  Moreover,  we  let  the  generic  7*  be  equal  to 


the  ‘critical’  damping  defined  as  follows 


7*  =  2  V m i  kl 


(15) 


The  equations  of  motion  (14)  are  associated  with  the  initial  conditions 


u 


N 


(t  =  to)  =  (ulN){k  1},  i  =  0, . . . ,  N  -  1;  u%(t  =  t0)  =  5 


ulN(t  =  t0)  =  0,  i  =  0, 


,N 


(16) 


where  (ulN)(k  ^  (i  =  1 , . . . ,  N)  are  the  displacements  of  the  masses  at  the 
external  time  r  =  iy_i .  Equations  (14),  (16)  are  numerically  integrated 


through  a  fourth-order  Runge-Kutta  integration  scheme,  up  to  an  internal 


time  t]  such  that  it  results  \al+l  —  al\  <  10  6|aA  |  for  all  i  E  {1, . . . ,  N  —  1}. 
When  the  internal  equilibrium  is  reached,  we  set  k  <—  k  +  1  and  re-iterate 


problem  (14). 


For  the  micro-meso  convergence  study  of  Section  5T  we  consider  the 
stress-strain  law  described  by  Equation  (J6])  and  illustrated  in  Fig.  [3j  for  each 
microscopic  spring.  For  the  simulation  of  compression  tests  on  real  CNT 
foams,  we  instead  adopt  nonuniform  chains  of  mesoscopic  spring  character¬ 
ized  by  a  suitable  numerical  regularization  of  the  stress-strain  pattern  shown 
in  Fig.  [4}  In  detail,  we  introduce  a  ‘hardening’  type  regularization  consisting 
of  the  following  stress-strain  law  (Fig.  [5]) 


(  a{a,i)  _ 
a  ~  (1 


<T  — 


a 


m  =  4  +  4+ (4 -4), 


(e®  <  £la)  and  (flag^  ^  7^  c), 

(£a  <  £*  <  £*)  and  (flag^_1)  =  a), 


crh’d  =  cr®  +  Act®  +  klh_(£z —  e^),  (ela  <  £®  <  el)  and  (flag(fc  =  c), 


(c,i)  =  kj^-ep 


{tc  <  £*)  and  (flag*-*’  ^  7^  a), 


(17) 
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where,  for  each  r  =  rk  (k  —  1, . . . ,  M),  we  set 


flag(fc) 


a ,  (e*  <  ejj  and  (flag^-1^  7^  c), 

c,  (ej.  <  e*)  and  (flag(fc-1)  7^  a), 

flag^-1-*,  ((£la  <  £*  <  £*)  and  (flag^1)  =  a)), 

or  ((ila  <  £l  <  il)  and  (flag('fc_1)  =  c)). 


(18) 


The  quantities  kl0,  Aa\  klc ,  e“,  £*,  klh+  and  k\_  in  (17)  are  constitutive 
parameters  (seven  independent  parameters),  while  £la  and  £*  are  computed 
by  solving  for  £l  the  equations 


=  cr(e’*) ,  a(c>i)  =  cx(cM)  (19) 

respectively.  One  gets  back  to  a  perfectly  plastic  mesoscopic  response  by 
setting  the  regularization  (‘hardening’)  parameters  klh+  and  klh_  equal  to 
zero.  The  quantities  £l  are  given  by  ([2])  under  the  replacement  of  u) v  with 
u%. 


O1 


Figure  5:  (Color  online)  Hardening-type  regularization  of  the  mesoscopic  response. 
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5.  Applications 


The  present  Section  deals  with  applications  of  the  above  numerical  model 
to  a  convergence  study  and  the  fitting  of  laboratory  tests  on  the  cyclic  com¬ 
pression  of  CNT  foams.  In  the  first  case,  series  of  springs  with  constitutive 
equations  ([6])  were  considered  (model  1:  micro- rneso  transition),  letting 
the  total  number  of  springs  of  the  system  increase  progressively  (micro-meso 


transition).  In  the  second  case,  spring  chains  with  constitutive  equations  ( 17) 


were  instead  examined  (model  2:  meso-macro  transition).  The  optimal  fit 
of  available  experimental  stress-strain  curves  was  performed  by  playing  with 
the  constitutive  parameters  of  equations  0  and  the  number  of  mesoscopic 
springs  N.  Depending  on  the  value  of  N ,  a  significant  number  of  parameters 
may  need  to  be  identified.  This  lead  us  to  employ  a  fitting  procedure  based 
on  genetic  algorithms,  which  have  been  proved  to  be  well-suited  for  multi¬ 


modal  non-convex  optimization  (cf.  e.g.  El  Sayed  et  al. ,  2008  Fraternali  et 


ah]  2009).  In  all  the  examined  examples,  we  name  global  strain  the  quantity 
e  =  (L  —  £)/L,  t  denoting  the  total  deformed  length  of  the  chain. 

5.1.  Convergence  study  and  micromechanics  of  ‘plastic’  steps 

We  examined  uniform  chains  with  fixed  length  L  =  860  fin i  and  increasing 
number  of  microscopic  springs  N,  subject  to  a  complete  loading- unloading 
compression  cycle  up  to  a  global  strain  e  =  0.85  (as  in  the  experiments  ex¬ 
amined  in  Section  5.2.2).  We  employed  model  ffl  with  kc  =  ko  (‘symmetric’ 


case);  e0  =  0;  and  the  material  properties  listed  in  Tab.  1  of  the  Appendix  for 
all  the  springs.  Fig.[6]shows  the  numerically  computed  dissipation  DN  versus 
strain  £  and  stress  a  versus  e  responses  of  such  a  model  for  different  numbers 
of  springs  ( N  =  5,  10,  15  and  50).  The  a-e  plot  in  Fig.  [6]  highlights  that 
the  global  stress-strain  response  alternates  elastic  steps  and  ‘plastic’  jumps 
of  a  at  constant  e,  oscillating  converging  to  the  ‘perfectly  plastic’  mesoscopic 
behavior  shown  in  Fig.  [4]  The  D^-e  plot  in  Fig.  [6]  instead  shows  that  D n 
converges  to  the  limit  dissipation  defined  in  equation  (13).  The  above  results 


therefore  confirms  the  theoretical  predictions  of  Section  [3]  Selected  equilib¬ 
rium  configurations  and  a  deformation  animation  (online  version  only)  of  the 
model  with  N  =  5  springs  are  given  in  Fig.  [7]  One  can  easily  recognize  the 
configurations  corresponding  to  the  plastic  steps  of  the  microscopic  chain, 
displaying  the  snap  of  a  single  spring  and  the  simultaneous  elastic  rearrange¬ 
ment  of  the  remaining  springs  (cf.,  e.g.,  configurations  3  and  5  from  the  left). 
It  is  worth  noting  that  the  ordering  of  the  springs  is  of  no  relevance  in  the 
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present  case.  Therefore,  the  succession  of  the  spring  snaps  can  indifferently 
proceed  from  the  bottom  to  the  top  of  the  chain,  as  shown  in  Fig.  [TJ.  from 
top  to  bottom,  or  in  random  sequence.  In  a  real-world  chain  of  springs,  the 
parameters  as  k0  or  ea  would  never  be  perfectly  identical  and  the  small  devi¬ 
ations  of  the  different  springs  would  break  the  symmetry  and  determine  the 
ordering  when  the  springs  snap. 


5.2.  Fitting  of  experimental  results  on  compressed  CNT  foams 

We  examined  the  experimental  results  of  a  cyclic  compression  test  on  a 
doubly  anchored  CNT  foam  run  at  the  Graduate  Aerospace  Laboratories  of 


the  California  Institute  of  Technology,  and  the  results  given  in  Cao  et  al. 


(2005)  on  the  cyclic  compression  of  a  foam-like  CNT  him.  We  fit  model 
ff  2  both  to  he  first  loading/unloading  cycle  of  the  examined  foams,  and 
to  selected  cycles  following  an  initial  mechanical  preconditioning,  employing 
the  Breeder  Genetic  Algorithm  (BGA)  presented  in  the  Appendix.  For  the 
fit  of  the  first  cycles  we  assumed  el0  =  0  in  each  spring,  while  for  the  cycles 
following  the  foam  preconditioning  we  set  el0  =  e o  in  each  spring,  eo  being  the 
macroscopic  (permanent  or  transient)  strain  observed  at  zero  stress  at  the  end 
of  the  previous  loading  cycle.  For  the  cycles  following  the  preconditioning, 
we  considered  two  different  fitting  models:  one  accounting  for  hysteresis 
(numerical  (1):  Act*  <  0  in  each  mesoscopic  spring),  and  the  other  without 
hysteresis  (numerical  (2):  A  a1  =  0  in  each  mesoscopic  spring).  Since  all  the 
adopted  fitting  models  account  for  macroscopic  time-independent  behavior 
(cf.  Sect.  [3]),  our  fitting  results  are  independent  of  the  strain  rate  actually 
applied  in  the  experiments. 


5.2.1.  Compression  tests  on  a  doubly  anchored  CNT  foam 

We  performed  cyclic  compression  tests  on  vertically  aligned  multi-walled 
carbon  nanotube  forests  (800 /im  in  length  with  sample  area  ~  14mm2)  grown 
by  chemical  vapor  deposition  (CVD)  using  ferrocene  and  toluene  as  precur¬ 
sors.  The  average  diameter  of  the  as-grown  CNTs  was  ~  50nm.  The  as-grown 
CNTs  foams  were  partially  anchored  between  two  polymer  layers  to  provide 
structural  support  and  sample  transportability. 

For  anchoring  the  CNT  foams  on  a  substrate  we  spin-coated  polydimethyl- 
siloxane  (PDMS)  on  top  of  a  glass  slide  at  800  RPM  to  get  50  —  100 pm  thick 
films.  Before  curing  the  polymer  we  partially  embedded  the  CNT-foams 
at  80°C  for  1  hour.  After  curing,  the  carbon  nanotubes  protruding  from 
the  substrate  showed  excellent  vertical  alignment  with  an  average  height  of 
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~  750/xm.  The  process  was  repeated  turning  the  sample  upside  down  to 
obtain  a  ‘doubly  anchored’  system  (i.e.  sandwich  structure  with  polymer  on 
both  sides  and  the  CNT  foam  in  the  middle).  Typical  experimental  results 
obtained  from  cyclic  compression  tests  are  shown  in  Fig.  I  (top). 

We  separately  fit  the  first  and  the  fourth  stress-strain  cycles  shown  in 
Fig.  [8]  to  model  ff  2,  considering  the  general  ‘asymmetric’  case  with  klc  7^  kl0 
(i  —  0, . . . ,  N  —  1)  and  a  chain  with  iV  =  4  springs.  For  the  fourth  cycle  we 
accounted  for  an  initial  strain  £0  =  0.20,  which  approximatively  corresponds 
to  the  strain  measured  at  zero  stress  at  the  end  of  the  third  loading  cycle 
(Fig.  [8j  top).  The  best  fit  parameters  obtained  through  BGA  optimization 
are  given  in  Tab.  2  of  the  Appendix.  A  comparison  between  best-fit  and  ex¬ 
perimental  overall  stress-strain  curves  is  shown  in  Fig. [8] (center  and  bottom). 
One  observes  that  the  non-uniform  dissipative  mass-spring  model  is  able  to 
capture  the  real  hysteretical  behavior  of  the  examined  CNT  foam,  both  at 
their  pristine  state  (first  cycle)  and  after  mechanical  preconditioning  (fourth 
cycle),  through  a  multi-plateaux  overall  stress-strain  profile.  In  particular, 
Fig.  (8}bottom  and  Tab.  2  of  the  Appendix  show  that  the  response  after 
preconditioning  can  be  roughly  described  through  a  suitable  non- dissipative 
model. 


5.2.2.  Compression  tests  on  a  foamlike  CNT  film 

We  fit  a  cyclic  compression  test  given  in  Cao  et  al.  (2005)  for  a  CNT  foam¬ 
like  film  with  thickness  L  =  860  /irn.  The  analyzed  experiment  performed 
1000  loading/unloading  cycles  up  to  a  global  strain  £  =  0.85  (cf.  Fig.  3B  of 
Cao  et  al. ,  2005'  and  Fig.  3  of  the  Appendix).  The  fitting  of  the  ‘symmet¬ 
ric’  formulation  of  model  fi  2  (klc  =  kl0  in  each  spring)  to  the  first  cycle  is 
presented  in  the  Appendix.  A  dramatic  improvement  of  the  fitting  ability 
of  model  2  was  obtained  by  considering  the  ‘asymmetric’  case  (/c®  7^  kfi). 
As  shown  in  Tab.  4  of  the  Appendix,  we  were  indeed  able  to  reduce  the 
fitting  fitness  to  0.808  MPa  by  using  a  BGA-optimized  ’’asymmetric”  mul¬ 
tiscale  model  with  10  springs.  Such  a  fitness  value  is  markedly  lower  than 
that  obtained  for  the  corresponding  ’’symmetric”  case  (1.462  MPa,  cf.  Tab. 
3  of  the  the  Appendix).  We  were  able  to  optimally  fit  also  the  1000th  cycle 
of  the  examined  experiment  to  the  ‘asymmetric’  model,  obtaining  a  fitting 
fitness  of  0.444  MPa  through  the  BGA-optimized  5  spring  model  described  in 
Tab.  4  of  the  Appendix.  The  excellent  match  between  numerical  and  exper¬ 
imental  results  for  the  present  case  is  illustrated  by  Fig.  [9j  which  compares 
predicted  and  measured  overall  stress-strain  responses.  For  the  first  cycle 
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we  observe  that  the  numerical  multi-plateaux  response  closely  approaches 
the  continuous  experimental  recording  as  the  number  of  mesoscopic  springs 
increases  from  5  to  10  (Fig.  [9j  top  and  center).  Regarding  the  1000th  cycle, 
we  assumed  e0  =  0.14  as  suggested  in  |Cao  et  al.  (2005).  Fig.  [9j  bottom 


shows  that  the  experimental  stress-strain  profile  corresponding  to  the  1000th 
cycle  is  already  excellently  approximated  by  a  5  spring  model  with  dissipa¬ 
tion  (numerical  (1)),  and  roughly  described  by  an  analogous  model  without 
dissipation  (numerical  (2)).  One  observes  from  Tabs.  3  and  4  of  the  Ap¬ 
pendix  that  the  average  value  of  the  stiffness  kl0,  among  all  the  springs,  is 


approximative^  equal  to  the  elastic  modulus  estimated  by  Cao  et  al.  (2005) 
for  the  present  CNT  foam  (^  50  MPa). 


6.  Discussion  and  Outlook 

In  this  article  we  derived  a  mechanical  model  describing  the  behavior 
of  CNT  foams  under  uniform  compression,  through  the  concept  of  bi-stable 


springs  presented  in  Puglisi  and  Truskinovsky  (2000,  2002,  2005),  which  has 


been  recently  applied  to  open-cell  polyurethane  foams  (Pampolini  and  Del 


Piero,  2008). 


The  proposed  model  differs  from  other  bi-stable  mass-spring  models  avail- 


able  in  the  literature  (refer  e.g.  to 

31esgen 

2007 

Br aides  and  Cicalese 

2007 

Charlotte  and  Truskinovsky,  2002 

2008f  Ortiz,  |1999|  Puglisi  and  r 

Tuski- 

novsky,  2000,  2002,  2005,  and  references  therein),  due  to  the  presence  of 
an  intermediate,  mesoscopic  scale,  placed  in  between  the  microscopic  scale 
of  the  bi-stable  springs,  and  the  macroscopic  scale  of  the  entire  structure. 
The  microscopic  scale  aims  to  describe  the  dynamic  snapping  of  the  carbon 


nanotubes,  due  to  local  buckling  (Cao  et  al.,  2005  Zbib  et  al.  2008  Misra 


et  al.  2009;  Hutchens  et  al.  2010),  while  the  mesoscopic  scale  is  intended 


to  describe  the  time-independent  hysteretic  behavior  of  finite  portions  of  the 
CNT  foam.  The  latter  typically  follows  from  kinking,  entanglement  of  the 
tubes,  friction  and  other  microscopic  dissipative  effects  (Suhr  et  al.  2007; 


Teo  et  al. 

2007; 

Tao  et  al. 

2008 

Deck  et  al. 

2007 

Misra  et  al. 

2009 

Another  relevant  difference  between  the  present  model  and  most  of  the 
available  bi-stable  mass-spring  models  consists  of  the  fact  that  the  current 
model  accounts  for  nonuniform  mesoscopic  spring  properties,  while  other 
models  usually  consider  uniform  chains  of  bi-stable  springs  (see  e.g.  Puglisi 


and  Truskinovsky  2000  2002;  Pampolini  and  Del  Piero ,  [2008 ) .  Such  a  me¬ 
chanical  inhomogeneity  allows  us  to  account  for  hardening  of  the  macroscopic 
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response  in  the  post-buckling  range  (Cao  et  al.,  2005  Misra  et  al.,  2009). 

The  numerical  results  given  in  Section  [5]  demonstrated  that  our  methods 
are  capable  of  recovering  available  experimental  results  on  the  compression  of 
CNT  foams  with  excellent  agreement.  We  fitted  the  parameters  of  the  meso¬ 
scopic  springs  through  a  Breeder  Genetic  Algorithm,  which  is  well  suited 
for  global  non-convex  optimization.  The  fitting  of  experimental  stress-strain 
curves  allowed  us  to  recognize  that  non-uniform  dissipative  mass-spring  mod¬ 
els  are  well  suited  to  capture  the  main  features  of  the  real  compressive  re¬ 
sponse  of  CNT  foams,  and  specifically  strain  localization  due  to  CNT  kinking 
and  time-independent  hysteresis.  The  latter  was  found  to  be  relevant  during 
loading/unloading  from  the  pristine  state,  and  progressively  decaying  after 
mechanical  preconditioning.  By  setting  to  zero  the  dissipation  of  the  model, 
we  were  led  to  obtain  a  rough,  non-linearly  elastic  approximation  of  the  ex¬ 
amined  experimental  behaviors  after  preconditioning. 

The  distinctive  features  of  the  multiscale  model  presented  in  this  work 
highlight  its  potential  use  to  describe  the  mechanical  behavior  of  multilay¬ 
ered  systems  composed  of  alternating  CNT  foams  and  anchoring  polymeric 
films  (Misra  et  ah,  2009).  We  will  address  such  an  extension  of  the  current 


model  in  future  work.  Other  future  directions  of  the  present  research  might 
regard  the  description  of  permanent  deformation  through  fatigue  or  damage 
mechanisms;  the  inclusion  of  rate-dependent  dissipative  effects  at  the  macro¬ 
scopic  scale;  and  the  modeling  of  long-range  interactions  mimicking  van  der 
Waals  forces. 
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Appendix.  Supplementary  data 


Supplementary  data  associated  with  Section  [5]  of  this  article  can  be  found 
in  the  online  version. 
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;ure  6:  (Color  online)  Response  of  uniform  chains  of  microscopic  bi-stable  springs  with 
iperties  shown  in  Tab.  1  of  the  Appendix:  D ^  dissipation  (kJ/m2);  a  global  stress 
Pa);  e  global  strain  (filled  marks:  loading;  unfilled  marks:  unloading). 
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tion) 


Figure  7:  (Color  online)  Selected  equilibrium  configurations  (left)  and  deformation  ani¬ 
mation  (right  -  online  version  only)  of  the  model  described  in  Tab.  1  of  the  Appendix,  for 
TV  =  5. 


a[|VPa] 
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Figure  8:  (Color  online)  Fitting  of  compression  tests  on  a  doubly  anchored  CNT  foam 
(top)  to  non-uniform  ‘asymmetric’  spring  models  (properties  in  Tab.  2  of  the  Appendix). 
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Figure  9:  (Color  online)  Fitting  of  experimental  compressive  stress-strain  curves  of  CNT 
foams  (Cao  et  al.  2005)  to  non-uniform  ‘asymmetric’  spring  models  (properties  in  Tab.  4 
of  the  Appendix). 
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Addendum  to  Section  5.  Applications 

Section  5.1.  Convergence  study  and  micromechanics  of  ‘plastic’  steps 

The  material  properties  used  for  the  convergence  study  presented  in  Sec¬ 
tion  5.1  of  the  paper  are  given  in  Tab.  [l] . 


k0  [MPa] 

Aa/aa 

£a 

£c 

50.00 

-0.7 

0.25 

0.77 

Table  1:  Mechanical  properties  of  uniform  chains  of  microscopic  bi-stable  springs  employed 
in  Section  5.1  of  the  paper. 


Section  5.2.  Fitting  of  experimental  results  on  compressed  CNT  foams 

The  fits  presented  in  Section  5.2  of  the  paper  consider  single  loading- 
unloading  cycles,  described  by  data  sets  of  the  form 


£r  £();  ^7/ r= 


i ,-,Nd  r  > 


(1) 


where  er  are  experimental  observations  of  the  global  strain  e  (hard-device 
conditions);  £q  is  the  permanent  strain  eventually  accumulated  during  a 
previous  load  history  (mechanical  preconditioning);  ar  are  the  experimen¬ 
tal  recordings  of  the  overall  stress  a;  and  Nj  is  the  number  of  data  points. 
We  look  for  the  best-fit  values  of  the  constitutive  parameters  of  model  ff  2 


P=  {{fco>  A<x\  K ,  4  kh+  4-L=i,...,tv-i}  , 

under  simple  bounds  of  the  form 

peD  =  \fi,p'*\x...x\$,pf\, 


(2) 

(3) 


where  P  is  the  overall  number  of  parameters.  We  set  P  =  51V,  by  prescribing 
the  ratios  k'h+/k0  and  klh_/k0.  In  detail,  we  fixed  klh+  =  k0  x  10“3  and 
klh_  =  k0  x  0.5  x  10~2  for  the  first  cycle  of  the  experiment  analyzed  in 
Section  5.2.1;  klh+  =  klh_  =  ko  x  10-5  for  the  fourth  cycle  of  the  experiment 
analyzed  in  the  same  Section;  and  kf  .  =  klh_  =  ko  x  10”2  for  the  experiments 
analyzed  in  Section  5.2.2  The  fitting  performance  of  a  given  set  of  parameters 
p  was  evaluated  through  the  fitting  fitness  function 


fip)  =  max  |  (Tr{p)  -  ar 

r=l,...,Nd 


(4) 


2 


which  is  the  maximum-norm  of  the  piecewise  continuous  residuals  a— a.  Here, 
a r(p)  denotes  the  numerically  predicted  overall  stress  for  e  =  er,  coinciding 
at  equilibrium  with  the  stress  in  each  individual  spring.  The  multivariate 
minimization  problem 

min  f(p),  (5) 

peD 


is  expected  to  be  strongly  non-convex  (Ogden  et  al. 

2004) 

and  well 

suited 

ployed 

for  genetic  algorithms  (Schmitt  2004  El  Sayed  et  al. 

2008 

).  We  em 

the  Breeder  Genetic  Algorithm  (BGA)  presented  in 

De  Falco  et  al. 

(1996) 

and  successfully  used  as  a  parameter  identification  tool  in 

Fraternal 

et  al. 

(2009).  We  used  a  population  size  of  2 P  individuals;  an  initial,  randomly- 
chosen,  truncation  rate  equal  to  15%,  extended  intermediate  recombination, 
mutation  rate  in  the  interval  [10%,50%],  and  a  maximum  number  of  gener¬ 
ations  equal  to  200.  We  refer  the  reader  to  De  Falco  et  al.  (|1996)  for  further 
technical  details  of  the  employed  BGA. 


Section  5.2.1.  Compression  tests  on  a  doubly  anchored  CNT  foam 

Tab.  [2]  illustrates  the  best  fit  material  parameters  obtained  for  the  exper¬ 
iment  illustrated  in  Section  5.2.1  of  the  paper. 


Section  5.2.2.  Compression  tests  on  a  foamlike  CNT  film 

We  fit  the  ‘symmetric’  formulation  of  model  ff  2  (klc  =  kf  in  each  spring) 
to  the  first  cycle  of  the  cyclic  compression  test  given  in  Cao  et  al.  (2005]), 
employing  chains  with  3,  5  and  10  mesoscopic  springs.  The  best  fit  pa¬ 
rameters  obtained  for  this  case  are  given  in  Tab.  [3]  and  the  corresponding 
stress-strain  plots  are  shown  in  Fig.  [l]  One  observes  that  the  matching  be¬ 
tween  predictions  and  experimental  recordings  appreciably  increases  by  pro¬ 
gressively  adding  mesoscopic  springs  to  the  model  under  consideration  (the 
fitting  fitness  function  /  decreases  from  1.586  MPa  to  1.462  MPa,  by  letting 
the  number  of  springs  increase  from  3  to  10).  A  good  matching  between 
theory  and  experiments  was  also  observed  for  what  concerns  the  localization 
of  the  CNT  deformation  during  the  buckling  (‘plastic’)  phase.  Fig.  [2]  shows 
selected  equilibrium  configurations  and  a  deformation  animation  of  the  best 
fit  model  with  5  springs.  It  is  worth  noting  that  the  succession  of  the  spring 
snaps  depicted  in  Fig.  [2]  qualitatively  reproduces  the  progressive  kinking  of 
the  tubes  observed  during  the  test,  which  is  clearly  described  by  Figs.  1,  2 
and  4A  of  Cao  et  al.  (2005).  The  best  fit  parameters  obtained  for  ‘asym¬ 


metric’  formulation  of  model  ff  2,  which  is  presented  in  Section  5.2.2  of  the 


3 


4  springs, 

1.  cycle ,  e0  =  0,  /  = 

:  0.491  MPa 

spring  # 

&o  [MPa] 

£i 

Gfl 

°c 

hi  /  hi 
"'c/  ^0 

3 

11.82 

-0.82 

0.09 

0.95 

9.06 

2 

48.89 

-0.88 

0.05 

0.93 

8.66 

1 

86.95 

-0.75 

0.06 

0.61 

9.58 

0 

16.92 

-0.79 

0.27 

0.56 

9.77 

4  springs, 

q/c/e,  £0  =  0.20,  numerical  (1),  /  =  0.484  MPa 

spring  # 

kl0  [MPa] 

ACT*/< 

ei 

°a 

< 

K/K 

3 

35.93 

-0.69 

0.03 

0.84 

92.73 

2 

10.76 

-0.98 

0.88 

0.88 

59.37 

1 

20.42 

-0.86 

0.19 

0.27 

25.93 

0 

7.16 

-0.16 

0.02 

0.95 

96.70 

4  springs, 

4-  cycle,  £0  =  0.20,  numerical  (2),  /  =  0.746  MPa 

spring  # 

kl0  [MPa] 

Aff’/ffi 

ei 

< 

K/K 

3 

76.29 

0.00 

0.62 

0.62 

84.16 

2 

10.76 

0.00 

0.88 

0.88 

59.37 

1 

1.66 

0.00 

0.73 

0.73 

76.13 

0 

2.24 

0.00 

0.13 

0.72 

92.55 

Table  2:  Mechanical  properties  of  ‘asymmetric’  spring  models  fitting  compression  tests  on 
a  doubly  anchored  CNT  foam  (cf.  Section  5.2.1  of  the  paper). 
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paper,  are  given  in  Tab.  |4j 


3  springs,  1.  cycle ,  /  =  1.568  MPa 


spring  # 

kl0  [MPa] 

A,t'K 

£i 

ca 

°C 

2 

52.30 

-0.78 

0.17 

0.78 

1 

52.30 

-0.78 

0.23 

0.88 

0 

52.30 

-0.78 

0.27 

0.73 

5  springs,  1.  cycle ,  /  =  1.512  MPa 


spring  # 

kl0  [MPa] 

Aa‘/a‘a 

£* 

°a 

£i 

°c 

4 

45.97 

-0.92 

0.24 

0.78 

3 

44.55 

-0.76 

0.23 

0.72 

2 

44.58 

-0.78 

0.21 

0.80 

1 

44.45 

-0.76 

0.30 

0.84 

0 

47.76 

-0.73 

0.25 

0.83 

10  springs,  1.  cycle ,  /  =  1.462  MPa 


spring 

kl0  [MPa] 

£* 

ca 

ei 

°c 

9 

46.76 

-0.81 

0.20 

0.86 

8 

58.04 

-0.86 

0.58 

0.54 

7 

51.92 

-0.80 

0.25 

0.87 

6 

48.27 

-0.88 

0.22 

0.87 

5 

60.76 

-0.86 

0.21 

0.60 

4 

58.10 

-0.72 

0.25 

0.81 

3 

53.87 

-0.70 

0.18 

0.80 

2 

44.27 

-0.66 

0.25 

0.89 

1 

59.50 

-0.64 

0.19 

0.64 

0 

59.73 

-0.76 

0.21 

0.88 

Table  3:  Mechanical  properties  of  ‘symmetric’  spring  models  fitting  results  by  Cao  et  al.| 
(2005). 
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5  springs,  1.  cycle ,  £0  =  0,  /  =  0.870  MPa 


spring  # 

kl0  [MPa] 

A  o'/K 

£i 

ei 

°c 

K/K 

4 

39.06 

-0.69 

0.28 

0.86 

4.70 

3 

43.41 

-0.65 

0.31 

0.81 

4.71 

2 

48.48 

-0.90 

0.21 

0.81 

3.90 

1 

38.15 

-0.73 

0.24 

0.76 

5.00 

0 

46.89 

-0.77 

0.26 

0.80 

4.33 

10  springs,  1.  cycle ,  e0  =  0,  f 

=  0.808  MPa 

spring  # 

kl0  [MPa] 

Aff’K 

el 

ca 

Cc 

hl  /  lcl 

9 

45.15 

-0.70 

0.24 

0.85 

1.85 

8 

48.16 

-0.56 

0.30 

0.66 

5.82 

7 

39.88 

-0.60 

0.25 

0.72 

7.91 

6 

41.73 

-0.61 

0.19 

0.76 

2.70 

5 

54.97 

-0.69 

0.21 

0.86 

7.58 

4 

59.29 

-0.92 

0.21 

0.86 

2.20 

3 

44.24 

-0.89 

0.24 

0.86 

4.97 

2 

48.96 

-0.79 

0.27 

0.88 

9.68 

1 

61.24 

-0.69 

0.16 

0.83 

7.16 

0 

43.92 

-0.67 

0.29 

0.83 

6.69 

5  springs, 

1000.  cycle ,  £q  —  07 

A,  numerical 

1),  /  =  0.444  MPa 

spring  # 

kl0  [MPa] 

A  <t‘K 

el 

Ga 

°c 

K/K 

4 

8.66 

-0.28 

0.19 

0.90 

27.35 

3 

2.55 

-0.27 

0.74 

0.74 

47.35 

2 

86.24 

-0.96 

0.59 

0.63 

49.84 

1 

57.13 

-0.31 

0.44 

0.60 

20.15 

0 

2.95 

-0.71 

0.49 

0.85 

38.13 

5  springs, 

1000.  cycle ,  £0  =  07 

A ,  numerical 

2),  /  =  0.926  MPa 

spring  # 

klQ  [MPa] 

Act’K 

°c 

K/K 

4 

4.40 

0.00 

0.37 

0.89 

45.19 

3 

54.24 

0.00 

0.44 

0.71 

35.32 

2 

4.62 

0.00 

0.59 

0.63 

26.10 

1 

1.83 

0.00 

0.87 

0.87 

32.26 

0 

11.38 

0.00 

0.09 

0.75 

33.84 

Table  4:  Mechanical  properties  of  ‘asymmetric’  spring  models  fitting  results  by  Cao  et  ah 
(2005). 
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Figure  1:  (Fitting  of  the  first  cycle  of  a  compression  test  on  a  foamlike  CNT  film  to 
non-uniform  spring  models  (properties  in  Tab.[3|. 
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Figure  2:  Selected  equilibrium  configurations  (left)  and  deformation  animation  (right)  of 
the  model  with  five  springs  described  in  Tab.  [3j 
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